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We analyze the class of non-linear electrodynamics minimally coupled to gravitation supporting 
asymptotically flat non Schwarzschild-like elementary solutions. The Lagrangian densities governing 
the dynamics of these models in fiat space are defined and fully characterized as a subclass of 
the set of functions of the two standard field invariants, restricted by requirements of regularity, 
parity invariance and positivity of the energy, which are necessary conditions for the theories to be 
physically admissible. Such requirements allow for a complete characterization and classification 
of the geometrical structures of the elementary solutions for the corresponding gravity-coupled 
models. In particular, an immediate consequence of the requirement of positivity of the energy is 
the asymptotic flatness of gravitating elementary solutions for any admissible model. The present 
analysis, together with the (already published) one concerning the full class of admissible gravitating 
non-linear electrodynamics supporting asymptotically flat Schwarzschild-like elementary solutions, 
completes and exhausts the study of the gravitating point-like charge problem for this kind of 
models. 

PACS numbers: 03.50.De, 04.40.-b, 04.70.Bw, ll.10.Lm 
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I. INTRODUCTION 



The study of gravitating non-linear electrodynamics 
(G-NED) generalizing the Reissner-Nordstrom (RN) so- 
lution of the Einstein-Maxwell field equations has re- 
ceived a considerable attention during the last decades, 
mainly due to the finding that (abelian and non-abelian) 
Born-Infeld (BI) models |l| arise, together with gravi- 
tation, in the low-energy regime of string and D-Branc 
physics @- In this context, asymptotically flat, static, 
spherically symmetric black hole solutions for the BI 
theory minimally coupled to gravity were obtained in a 
number of papers [3|. Moreover, other G-NED models 
supporting electrically charged black hole solutions have 
been obtained. As examples let us mention: the grav- 
itating generalization of the Euler-Heisenberg effective 
Lagrangian of Quantum Electrodynamics the gravi- 
tating versions of a Bl-like family of Lagrangians [5[ and 
of the logarithmic Lagrangian of Ref. [6fl ; the black holes 
generated by Coulomb-like fields in (2+1) dimensions Q; 
those associated to models preserving the conformal in- 
variance of Maxwell theory in any dimension Q as well 
as to models whose Lagrangian densities are defined as 
powers of the Maxwell Lagrangian [9(. 

Most of the examples considered in this context in the 
literature are particular cases of the full class of admissi- 
ble non-linear electrodynamics (NED) in flat space whose 
Lagrangians are defined as arbitrary functions of the two 
first-order field invariants. Admissibility amounts to the 
restriction of these functions by some conditions endors- 
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ing the physical consistency of the associated models (es- 
sentially, the positive definite character of the energy and 
the parity invariance, aside from some requirements for 
the proper definition and regularity of the Lagrangian 
functions and their elementary solutions). The sub- 
class of these models supporting stable and finite-energy, 
electrostatic spherically symmetric (ESS) solutions (non- 
topological solitons) in flat space has been fully charac- 
terized and extensively analyzed in Refs. (loj and fill ]. 
The ESS solutions of the models belonging to this sub- 
class behave at the center (r = 0) and asymptotically 
(r — > oo) in such a wa y th at the integral of energy is con- 
vergent. Next, in Ref. (T^, we have extended the analysis 
to those admissible models supporting ESS solutions for 
which the integral of energy may diverge at the center. 
The main purpose in this reference was the study of the 
gravitating electrostatic spherically symmetric (G-ESS) 
solutions of the family of G-NEDs obtained from the min- 
imal coupling of this class of admissible NEDs to gravity. 
The geometrical structure of the corresponding G-ESS 
solutions was exhaustively characterized in terms of the 
central and asymptotic behaviours of the ESS solutions 
(or, equivalently, in terms of the behaviours of the La- 
grangian functions in vacuum and on the boundary of 
their domains of definition at vanishing magnetic field). 
Aside from naked singularities, extreme black holes or 
two-horizons black holes, as in the RN case, other grav- 
itational configurations, such as extreme or non-extreme 
black points and single-horizon black holes, appeared for 
the models with ESS soliton solutions in flat space. They 
are in all cases asymptotically flat configurations with 
well defined Arnowitt-Deser-Misner (ADM) masses. 

In order to exhaust the analysis of the geometrical 
structure of G-ESS solutions associated to the full class 
of admissible NEDs minimally coupled to gravity we 
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must study (besides those considered in Ref.[T2|) the 
models supporting ESS solutions whose integral of en- 
ergy diverges asymptotically. This is the main pur- 
pose of the present paper. As we shall see, these G- 
ESS solutions exhibit rich geometrical structures and, de- 
spite their unpleasant asymptotic behaviour in absence of 
gravity, they are always asymptotically flat, even though 
they approach flatness at large distances slower than the 
Schwarzschild field, and the ADM mass cannot be de- 
fined. We shall use the adjective "anomalous" in calling 
this behaviour. 

All the above gravitating solutions contain curvature 
singularities at the center, event though the metric may 
be finite there. This is a consequence of a non-existence 
theorem formulated for purely electrically charged solu- 
tions jl 31 ] . Strictly speaking this theorem concerns the G- 
ESS solutions of NEDs with Maxwellian weak-field limit. 
Nevertheless, in the last section, we shall extend this 
theorem by proving the unavoidable presence of a cur- 
vature singularity at the center of the G-ESS solutions 
associated to any admissible NED model minimally cou- 
pled to gravity. Such a theorem can be circumvented for 
other kinds of configurations, such as purely magnetically 
charged solutions 1141 or through a coupling between dif- 
ferent structures [151 ] . Let us mention that some models 
containing singularity-free electrically charged black hole 
solutions, found in Ref. |16| . correspond to non-admissible 
Lagrangian densities which suffer "branching" features as 
functions of their arguments, as explained in 17]. 

In section II we summarize the main results on ad- 
missible NEDs in flat space, extending the classification 
of Refs.pjJ and 12 1 in order to include the models sup- 
porting asymptotically energy-divergent ESS solutions. 
Section III is devoted to the outline of the G-NED prob- 
lem in the ESS case and its resolution for the different 
types of models in terms of the boundary and vacuum be- 
haviours of the NED Lagrangians in flat space, focusing 
on the detailed analysis of those leading to asymptotically 
anomalous G-ESS solutions. We conclude in section IV 
with a summary and some perspectives. 



II. NON-LINEAR ELECTRODYNAMICS IN 
FLAT SPACE 

Let us recall the main results concerning the families 
of NEDs analyzed in references [ll[ and [l2[ and extend 
them to the full class of admissible NEDs. The gener- 
alized Lagrangian densities giving the dynamics of the 
fields are assumed to be general functions ip(X, Y) of the 
field invariants 



in the usual way. These functions are restricted by some 
"admissibility" conditions (such as regularity, positivity 
of the energy, parity invariance, etc.), which are necessary 
in order to define physically meaningful theories (Tlj . 

In terms of (p(X,Y), the symmetric (gauge- invariant) 
energy-momentum tensor reads 



and the energy density takes the form 



P = T m = 2^E 2 + 2^E-H-p(X,Y). (3) 

The necessary and sufficient condition for this energy 
functional to be positive definite and vanishing in vac- 
uum, can be shown to be [nfl 



<p(X,Y)>0, (4) 



which requires <p(Q,0) = 0, §^\(x>o,y=o) > and 
<p(X > 0, 0) > 0. The Euler equations associated to a 
Lagrangian density p(X, Y) are 



dip 
dX 



dp 

w 



F* 



= 0, 



(5) 



and the search for ESS solutions of the form E{r) 
E{r)*-,H = lead to the first-integral 



2 E(r) 



d<p 
dX 



Y=0 



= Q, 



(6) 



where Q is an integration constant, identified as the 
source electric charge. Owing to the admissibility con- 
dition <px\y=o > the signs of E(r) and Q are the 
same and, consequently, we can consider only the case 
Q > without loss of generality [TlJ] . The solutions with 
Q < are straightforwardly obtained through the re- 
placement E(r,Q,fJ,) — ► sign(Q)E(r,\Q\, fi). The form 
of the ESS solution of charge Q obtained from Eq.© 
(with E(r, Q) = vX; see Eqs.Q) is completely specified 
once the explicit expression of <p(X, 0) is given. Owing to 
the admissibility conditions, E(r, Q) is a monotonic func- 
tion of r [ll|. From Eq.Q the energy density associated 
to these solutions reads 



X = ~F I1U F> 1V = E 2 - H 2 
Y = -If^F*^ =2E-H, 



(1) 



where the field strength tensor F^ v and its dual F*^ v , as 
well as the electric E and magnetic H fields are defined 



, dp 
' dX 



Y=0 



E 2 (r,Q)-p(E 2 (r,Q),Y = 0), (7) 



whose integral in space, if finite, gives the total energy 
as a function of the charge. Obviously, the finite or 
divergent-energy character of the ESS solutions for a 
given model (aside from its linear stability) is related 
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to the form of the Lagrangian density ip(X, Y) governing 
the dynamics. The complete characterization of the fam- 
ilies of admissible NEDs supporting ESS non-topological 
soliton solutions in flat space has been performed in Refs. 
[Tol Hlj |. Besides these models, we are interested here in 
the whole family of admissible NEDs. As we shall see, 
when minimally coupled to gravity, they lead to G-ESS 
solutions which are asymptotically flat space-times. Let 
us consider the behaviour of the ESS field around the 
center and as r — > oo, assumed to be of the form [23| 

%4 0)~r", (8) 

and 



three subclasses: Bl models with — 2 < q < — 1 (ESS 
solutions damped slower than the coulombian field); B2 
models with q = —2 (coulombian damping); B3 mod- 
els with q < —2 (faster than coulombian damping). For 
models whose ESS solutions are damped with an expo- 
nent — 1 < q < in Eq.©, the integral of energy di- 
verges asymptotically. We shall call them IRD models 
(acronym for "infrared divergent"). 

Let us analyze now the behaviour of the Lagrangian 
densities around the values of their arguments corre- 
sponding to the central and asymptotic regions of the 
ESS solutions (see Fig.l). From the first-integral © we 
see that the Lagrangian density behaves, in the Al and 
UVD cases, as 



E(r — !> oo) 



(9) 



Using Eqs.Q and dU) we see that the positivity of the 
energy requires p < and q < 0. When p < in ((SJ) the 
ESS fields diverge at r — 0, whereas for p = the field 
must behave around the center as 



E{r -> 0) - a - 6r CT , 



(10) 



where the maximum field strength (a) and the exponent 
(a > 0) are characteristic parameters of the model, while 
the coefficient b is related to the charge of the particular 
ESS solution (the quantity a = Qb 2 l a turns out to be a 
universal constant for a given model |l2|). Consequently, 
the solutions \E(r,Q)\ of the admissible models decrease 
monotonically from +oo (if p < 0), or from a finite value 
\a\ (Hp — 0), at r = 0, and vanish asymptotically as 
r — > oo. The behaviour of the energy density around the 
center is given by 



<p(X,Y = 0) - V 7 , 
around X — E 2 (r — > 0) — > oo, where 
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p-2 
2p ■ 



(13) 



(14) 



The condition (fj| for the positivity of the energy reads 
in these cases 



(27 - l)X-< > 0, 



(15) 



and this equation confirms that if p > (so that the 
exponent 7 < 1/2) the energy density becomes negative 
around the vacuum. The Al cases correspond to 7 > 3/2 
whereas the UVD cases correspond to 1/2 < 7 < 3/2. 
For p — (case A2) the field at the center behaves as 
in Eq. (|10j) and the Lagrangian density behaves around 
(X = E 2 (r = 0) = a 2 , Y = 0) as 



p(r -> 0) - r v 



(11) 



<p(X,Y = 0) 



2aa 



(a 



(16) 



(see Eqs.© and (J7J) so that the integral of energy con- 
verges there if — 1 < p < and diverges if p < — 1. 
Following our previous conventions [III we denote class- 
Al models those for which the ESS solutions behave 
as in Eq.© with the exponent — 1 < p < and class- 
A2 models those supporting ESS solutions which behave 
as in Eq. (fT0)) (corresponding to p = in Eq.flS)). In 
both cases the integral of energy converges at the center. 
On the other hand we shall call "ultraviolet divergent" 
(UVD) those models with p < —1, for which the integral 
of energy of the ESS solutions diverges at the center. 
Asymptotically, the energy density behaves as 



if u 7^ 2. Here a and A are universal constants of the 
model UM. If a = 2 we have 



(p(X,Y = 0) aln(a- 



(17) 



Consequently, for a < 2 the Lagrangians exhibit a ver- 
tical asymptote at this point. On the other hand they 
take a finite value there with infinite slope for a > 2 (the 
BI model belongs to this case with a = 4). 
The same equation (I13[) . now with 



p(r — > 00) ~ r q 2 , 



(12) 



and the integral of energy converges there if q < —1, 
but diverges for — 1 < q < 0. In our previous work the 
asymptotically energy-convergent cases were split into 



7 



2q ' 



(18) 



describes the behaviour of the Lagrangian density around 
the vacuum (X = E 2 (r — > 00) = 0, Y = 0) and the 
ranges of values of the exponent are now 1 < 7 < 3/2 
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p(X,Y=0) 





B3 (q<-2) 
B2 (q=-2) 
1 (-2<q<-l) 

IRD (-l<q<0 




<-l 



e(Q) = £m(oo, Q) = e ex (0, Q) = e in (r, Q) + e ex (r, Q), 

(21) 

where e(Q) is the total energy of the ESS solution of 
charge Q, scaling as e(Q) = Q 3 ^ 2 e(Q = 1). 

The full set of admissible non-linear electromagnetic 
field theories is now obtained by including, besides the 
soliton-supporting ones, the families of admissible models 
whose associated ESS solutions are energy-divergent, ex- 
hibiting UVD, IRD or both behaviors. This extended 
set is now naturally classified into twelve families, by 
combining the three central- field behaviors (Al, A2 and 
UVD) with the four asymptotic behaviors (Bl, B2, B3 
and IRD). 



FIG. 1: Behaviours of the Lagrangian densities of admissible 
NEDs around the vacuum (cases Bl, B2, B3 and IRD) and 
on the boundary of their domain of definition at Y = (cases 
Al, A2 and UVD). In between both domains, these functions 
must be strictly monotonic increasing in order for the energy 
to be positive definite. The parameters p, a and q correspond 
to the exponents of r in Eqs.©, (|10[) and ((9J, respectively, 
determining the central and asymptotic field behaviours of 
the ESS solutions. 



in case Bl, j = 1 in case B2, 1/2 < 7 < 1 in case B3 
and 7 > 3/2 in case IRD. All these behaviours of the 
admissible NEDs are plotted in figure 1. 

We see from this discussion and the regularity of 
E(r > 0, Q) that we can define, in cases Al and A2, 
the internal-energy function 



,(r,Q) -4vr / R 2 p{R,Q)dR, 
Jo 



(19) 



giving the energy of an ESS solution of charge Q inside a 
sphere of radius r. Owing to the admissibility conditions, 
this is a monotonically increasing and convex function of 
r. This function cannot be defined in the UVD cases. In 
the asymptotic cases Bl, B2 and B3 we can define the 
external-energy function 



e ex (r,Q)=4iT R 2 p(R,Q)dR, 



(20) 



which gives the energy of the field outside the sphere of 
radius r. This is a monotonically decreasing and concave 
function [l2|], which cannot be defined in the IRD cases. 

The six combinations of cases Al and A2 for the 
central-field ESS behaviour with cases Bl, B2 and B3 for 
the asymptotic behaviour, lead to six families of models 
which exhaust the class of admissible NEDs supporting 
finite-energy ESS solutions. In such theories ei n (r,Q) 
and e e x(^j Q) are well defined and satisfy 



III. GRAVITATING NON-LINEAR 
ELECTRODYNAMICS 

Let us consider now the minimal coupling of admissi- 
ble NEDs to the gravitational field, defined through the 
action 



S — Sg + Sned 



R 



16ttG 



<P(X,Y) 



,(22) 

where, as usually, g and R are the metric determinant 
and the curvature scalar, respectively. The variation of 
this action with respect to the electromagnetic fields lead 
to field equations which are obtained from Eqs.© by the 
replacement of partial derivatives by covariant deriva- 
tives. The expression of the mixed components of the 
energy-momentum tensor of the electromagnetic fields, 
from which the Einstein equations can be immediately 
written, remain the same as in flat space in the static 
spherically symmetric case, where these equations, to- 
gether with the relations between the components of the 
energy-momentum tensor 



T ° = Tl = 2E 



2 d( P 
dX 



Y=0 



To 



-ip, (23) 



allow for the introduction of an adapted coordinate sys- 
tem where the interval takes the Schwarzschild-like form 



ds 2 = \{r)dt 2 



dr 2 
A(r) 



r 2 dQ 2 , 



(24) 



The field equations for the 



with dfl 2 = d9 2 + sin 2 
G-ESS solutions in this coordinate system take the same 
form as those in spherical coordinates in the flat-space 
problem for the ESS solutions. Consequently, the first- 
integral ([5]) and the form of the ESS field solutions, as 
functions of the radial coordinates r, are the same in 
both cases. Thus, the analysis and classification of the 
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ESS solutions of the admissible models in flat space in 
terms of their central and asymptotic behaviours, can be 
immediately translated to the gravitating problem. 

On the other hand, the Einstein equations in this G- 
ESS case take the form 



— (rA(r) - r) 
|^A(r)) 



-87rr 2 T n ° 



-8nr 



1 1 2^E? 



dX 



(25) 



whose compatibility can be easily established. The gen- 
eral solution of these equations reads 



X(r,Q,C) = l + 



C 
r 



8tt 



r z T»{r,Q)dr, (26) 



where C is an arbitrary constant, absorbing the constant 
of the indefinite integral of the energy density term in 
this formula. When the asymptotic behaviour of the ESS 
solutions belongs to cases Bl, B2 or B3, Eq. (|2l))) can be 
written as 



X(r,Q,M) = l- 2 Jl + 2 -^l, (27) 
r r 

where e ex {r, Q) is the external-energy function in flat 
space, defined in Eq. ff20|) . and the constant M is iden- 
tified as the ADM mass. The structure of the G-ESS 
solutions of these admissible NEDs minimally coupled to 
gravity has been classified and extensively analyzed in 
Ref.pjj]. This analysis concerned the models obtained 
from the UVD, Al and A2 central field behaviours com- 
bined with the Bl, B2 and B3 asymptotic behaviours. 
The corresponding metrics are asymptotically flat, be- 
having for large r as the Schwarzschild one, no matter 
the central-field behaviour. If this central behaviour is 
UVD, the structure of the metrics (f2~T|) for different val- 
ues of M and Q is qualitatively the same as that of 
the Reissner- Nordstrom solution of the Einstein-Maxwell 
field equations (naked singularities, extreme black holes 
or two- horizons black holes) . If the central behaviour be- 
longs to cases Al or A2 we have, in addition, new kinds 
of solutions (non-extreme single-horizon black holes and 
extreme and non-extreme black points). However, as al- 
ready mentioned in the introduction, the analysis of ad- 
missible models supporting IRD solutions in flat space, 
which lead also to asymptotically flat G-ESS configura- 
tions, was not considered in that reference. We shall now 
tackle this issue by exploring the properties of such so- 
lutions. These admissible models support ESS solutions 
in flat space behaving asymptotically as in Eq.© with 
> q > — 1 and, as a consequence of the asymptotic be- 
haviour of the energy density (|12[) , the integral of energy 
diverges at large r. Consequently, e ex (r,Q) in (|2T)|) can- 
not be defined. We are thus lead to three new cases by 
combining this asymptotic behaviour with the Al, A2 or 




FIG. 2: Qualitative behaviour of the energy curves y = 
£in(r, Q) (for fixed Q) for the different admissible IRD models 
(the curve UVD corresponds to the primitive of r 2 To(r, Q), 
defined up to an arbitrary constant). The slopes of these 
curves at r — diverge in the Al and UVD cases and take fi- 
nite values (= 8nQa) in the A2 cases. In all cases these curves 
are monotonic, convex and diverge asymptotically with van- 
ishing slope. The cutting points between these curves and the 
beam of dashed lines y = define the horizons of the con- 
figurations. The tangency points correspond to extreme black 
hole and extreme black point configurations. Otherwise the 
configurations are naked singularities and one or two-horizons 
black holes. The cutting points between the energy curves 
and the straight line C = provide the radius of the critical 
configurations (which are meaningless for the UVD case). 



UVD central- field behaviours, for which the methods of 
Ref.[l2| cannot be immediately applied. Let us analyze 
these cases separately. 



A. IRD-A1 case 

The exponents in Eqs.© and ([9]) range now in the 
intervals — 1 < p < and — 1 < q < 0, respectively. The 
integration of the Einstein equations (|23|) leads to 



X(r,Q,C) = 1 + 



C 2e m (r,Q) 
r r 



(28) 



where Ei n {r,Q) is the interior integral of energy in flat 
space, defined by Eq. (|19p and C is an arbitrary integra- 
tion constant. Owing to equation (fL2|) . et n (r, Q) diverges 
asymptotically as 



e in (r oo, Q) — r 



9+1 



oo, 



(29) 



exhibiting a horizontal parabolic branch. Thus the last 
term in the r.h.s. of (|2"5f dominates over the second one 
for large r and the metric approaches asymptotic flatness 
as 1 — A(r — > oo) ~ r q , slower than the Schwarzschild 
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field. Consequently, there is not a well defined mass for 
these configurations (in fact, the expression for the ADM 
mass [HI diverges for all the IRD cases) . At r — equa- 
tions (fTTj) and ([TO]) show that £i„(r = 0, Q) = and 
this function exhibits a vertical slope there (see figure 2). 
The term C/r in (|28|) dominates at small r and the met- 
ric function diverges at the center with the sign of C (see 
figure 3). The condition X(r,Q,C) = 0, determining the 
radii of the horizons, leads to 



A(r) 



(r,Q) 



C 



(30) 



which relates these radii and the values of C for a given 
charge. Consequently, the horizon locations are given 
by the intersection points between the beam of straight 
lines y — and the monotonically increasing func- 

tion y — Ei n (r,Q). The straight line of the beam which 
is tangent to the internal-energy function defines an ex- 
treme black hole configuration, which is determined from 
Eg. (I3TJ1) . aside from the condition 



dr 



47rr 2 T °(r,Q) 



(31) 



This condition gives the tangency point location 
Thextr(Q), which is the horizon radius of the extreme 
black hole configuration and, in addition, Eq. (|30|) gives 
the value of the constant C ex tr(Q) of this tangent line, 
both parameters being functions of the charge. There are 
now several possibilities, depending on the values of C as 
compared to the one of C extr {Q) (see figures 2 and 3): 

1) C > C extr (Q)' There are no horizons and the solu- 
tions exhibit naked singularities at the center, where 
the metric function A(r, Q) diverges to +oo (line I on 
Fig.3). 

2) C = Cextr(Q)- The solutions are extreme black 
holes exhibiting a degenerate horizon at rh ex tr(Q) (line 
II on Fig.3). 

3) < C < C ex tr{Q)- The solutions are two- horizons 
black holes. As C decreases in this range, the inner 
(Cauchy) horizon radius decreases from rh ex tr{Q) to zero, 
whereas the outer (event) horizon radius increases from 
i"he X tr{Q) to a value rh cr it{Q) (lines III and IV on Fig.3). 
For the limit of these solutions obtained as C — > + the 
metric function \(r,Q,C — > + ) diverges at r = 0, step- 
ping from +oo to -— oo in crossing a vanishing-radius inner 
horizon. 

4) C = 0: Consistently with the conventions of Ref.[l2T| 
we shall call this one the critical configuration. The met- 
ric function coincides with that of the limit of the so- 
lutions of the precedent case 3) for r > 0, but diverges 
to — oo at r = (no step at the center). There is a 
unique horizon at r = rhcrit{Q) and the configuration is 
a single-horizon black hole (dashed line V in Fig.3). 

5) C < 0: The configurations are single-horizon 
black holes and A(r, Q, C) increases monotonically from 
—oo atr = 0toA— s-lasr— ^oo (line VI on Fig.3). The 




FIG. 3: Behaviour of A(r) for the IRD-A1 configurations. 
When C > there can be naked singularities (curve I, 
C > Cextr(Q)), extreme black holes (curve II, C = C e xtr{Q)) 
or two-horizons black holes (curves III and IV, < C < 
C e xtr(Q)). Solutions with C < lead to single-horizon black 
holes (curve V, C = 0, corresponding to the critical configu- 
ration and curve VI for C < 0). The structure of the limit 
configurations obtained as C — > ± can be easily visualized 
in this figure. All curves approach asymptotic flatness slower 
than the Schwarzschild solution. Curves I to IV describe also 
the typical qualitative behaviour of the metric function in the 
IRD-UVD cases. 



limit of these solutions obtained as C — >• 0~ coincides 
with the critical solution of the precedent case for any r 
and, for r > 0, with the limit solution of case 3) above. 

The causal nature of the central curvature singularities 
of the 1) to 3) configurations (C > 0) is timelike whereas 
for both 4) and 5) configurations (C < 0) these singu- 
larities are spacelike. As mentioned above, the metric 
function of the limit configuration obtained as C — > + , 
coincides with the one of the critical configuration for 
r > 0, but differs by an infinite step at the center. This 
step introduces supplementary 5-like distribution terms 
in the curvature tensor which modify the causal structure 
of the singularity [l8[ . 

As a simple example illustrating these IRD-Al models 
let us look for a family of NEDs supporting ESS solutions 
of the form 



E(r,Q,s) 



W2 



(32) 



where < s < 1 parameterizes the family and Q is the in- 
tegration constant characterizing the solution, which will 
be identified as the electric charge. It is introduced under 
this form because the ESS solutions must depend on the 
radius through the ratio r / \J~Q, as results from the first- 
integral ((6]) (see Ref.[ll[). By using this first-integral we 
are easily lead to a family of NED Lagrangians <p(X, Y, s) 
behaving on the Y = axis as [24[ 
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<p(X,Y = 0,a) = Xi, 



(33) 



for X > 0, and where the parameter 7 = character- 
izing the different models in the family, is restricted to 
the range 3/2<7<oo. By integrating (J7J in this case, 
the internal-energy function takes the form 



t (r,Q,s) = 87r6(Q, S )r 



l-s 



(34) 



where 



Q(Q,s) 



2sQ_ 
s(l - s) V« + 2 



(l+s/2) 



(35) 



As expected this is a monotonically increasing and convex 
function of r, with a divergent slope at the origin and 
diverging with vanishing slope as r — > 00. The metric 
function reads 



A(r, C, Q, s) = 1 + — 
r 



16tt 



e(Q,«) 



(36) 



and defines a beam of curves having the form shown in 
Fig. 3. In all cases these metrics approach asymptotic 
flatness slower than Schwarzschild or RN configurations. 
The radii of the horizons and the values of the constants 
C(Q) for extreme black hole solutions are obtained from 
Eqs. (|30l) and (|31[) . respectively. The horizon radii of the 
critical black hole configurations result from the condi- 
tion A(r, C = 0, Q, s) = 0. These values are 



r hextr (Q,s) = (32tt(1 -s)Q(Q,s)) 1 / s 



Cextr{Q: s ) 



1 - S 



'hextr (Q,s) 



(37) 



t s\ \ Thextr 

rhcrit(Q,s) = — rrj— > r hextr (Q,s). 

(1 — s) L ' s 



The verification for this family of the remaining results 
l)-5) established for the general IRD-A1 models is now 
straightforward. 



B. IRD-A2 case 

In this case the expression of the metric function is still 
given by Eq. (|2"5]) . but £i n (r, Q) has now a finite slope at 
the center given by 



dei 



Or 



Airr 2 TS 



(38) 



as results from Eqs.© and (fTTj) (see Fig. 2). Conse- 
quently, the internal-energy function Si n (r,Q) takes a 



form which is similar to that of the precedent case IRD- 
Al, excepting that the slope at the center is now finite. 
In this case the metric function and its derivative behave 
near the center as 



C 

X(r -> 0,C,Q,a,b,a) ~ 1 - IBttQa H 



32irbQ 



<9A 
Or 



C 



r->0 r z 

for (7^2 and 



(a + l)(2~a) 
(<r + l)(2-a)'' 



C 



r a + Ar 2 (39) 



2Ar 



A(r 0,C,Q,a,b) - l-16nQa+- + 



^r 2 (l-21n(r))+Ar 



dX 
dr 



C 32nbQ . . 

— o rlll ( r ) 



2Ar, 



for a = 2, A being an integration constant. When a > 2 
this constant is given by A = ^j^-ip(a 2 , 0). Otherwise the 
value of A is not relevant because the associated term is 
not dominant in these equations. We must now distin- 
guish three subcases, according to 16nQa = 1 (see figure 
2). 

a) For the first subcase (167rQa < 1; curve A2b in 
Fig. 2) the beam of straight lines with C < cut the 
curve £in{r,Q) once, leading to single-horizon black 
hole configurations. When C > the configurations 
are timelike naked singularities for which the metric 
function diverges at the center if C > 0, and attains a 
finite positive value there (A(0) = 1 — IdirQa > 0) for 
the critical configuration (C = 0). For C < the con- 
figurations have a spacelike curvature singularity at the 
center, where the metric function is negatively divergent. 
The limit configuration (C — > CP) is a black point, for 
which the metric function steps from —00 at the center 
to A(r — > + ) — > 1 — WirQa > on the outer side of the 
vanishing-radius horizon, coinciding with the critical con- 
figuration for r > 0. The limit configuration (C + ) is 
a naked singularity, exhibiting also a singularity of the 
metric function A(r) at the center, which steps there from 
+00 to the finite value 1 — IQnQa. For r > this limit 
configuration coincides with the critical one and with the 
limit black point one. As easily seen from Eqs. (|3"9"]) and 

the slope of the critical metric function at the center 
is related to the value of the a parameter of the model, 
diverging for a < 1 (see Fig. 4), vanishing for a > 1 (see 
Fig. 5) and taking a finite value for a = 1 (see Fig. 6). 

b) For the second subcase (167r<5a > 1; curve A2a in 
Fig. 2) there is a tangency point of one of the straight lines 
of the beam ' ! -^- , defining extreme black hole config- 
urations, where C = C extr (Q) and the extreme horizon 
radius rh ex tr{Q) are given by equations (|30D and (|31l) . 
For different values of the constant C, corresponding to 
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the different straight lines of the above beam, the config- 
urations are similar to those of the IRD-Al case (C > 
Cextr(Q)'- naked singularities; < C < C ex t r (Q)- 
two-horizons black holes; C < 0: single-horizon 
black holes). The causal structure of the central singu- 
larity is timelike in the cases C > 0(A(r — > 0) — > +oo) 
and spacelike in the cases C < 0(A(r —> 0) — > — oo). For 
the critical configuration (C = 0) the metric function in 
Eqs. (|39f and (|40j) attains a negative finite value at the 
center, given by A(r = 0, Q) = 1 — WnQa < and, 
consequently, the causal structure of the central singu- 
larity is spacelike. For the limit configurations obtained 
as C — > ± the metric functions step at the center from 
±oo, respectively, to this value, coinciding with the crit- 
ical configuration for r > 0. Thus, the causal structure 
of the central singularities of these limit configurations is 
timelike and spacelike, respectively. As in the precedent 
subcase, the slope of A(r) at the center for the critical 
configuration depends on the value of the model param- 
eter a in the same way (see Figs. 4, 5 and 6). 

c) For the third subcase (WirQa = 1; curve A2c in 
Fig. 2) the straight line of the beam corresponding to 
C = is tangent to the internal-energy curve at r = 0. 
The corresponding critical configuration is an extreme 
black point. The metric function for this configuration 
vanishes at the center and its slope there depends on the 
value of the parameter a, as in the precedent subcases. 
The configurations with C > are naked singularities, 
whereas those with C < are single-horizon black 
holes. The limit configuration obtained as C — » + is 
a black point with a discontinuity of the metric func- 
tion, which steps from +oo at the center to zero on a 
vanishing-radius emerging horizon and coincides with the 
critical one for r > 0. The limit configuration obtained as 
C —> 0~ is a black point for which the metric function 
steps from — oo to zero in crossing the vanishing radius 
horizon, coinciding also for r > with the critical one. 
There is a null curvature singularity at the center of the 
extreme black point critical configuration, but the singu- 
larities of the black point limit configurations obtained 
as C — > ± are timelike and spacelike, respectively. 

As an illustrative example of this IRD-A2 case let us 
consider a one-parameter family of models supporting 
ESS solutions of the form 



A(r) 



£(r,Q, M ) = [(r 2 /Q) + M : 



21 -V2 



(41) 



where \i is the model parameter and Q > is the charge 
associated with the solution. Near the center these fields 
behave as 



E(r->0,Q,/x) -— , (42) 

and we see that the maximum field strength of Eq. (fTU|) 
is a = jtz -1 , the coefficient b = (2Qfi 3 )" 1 and a — 2. The 
use of the first-integral ^ leads to the expression 







C>0 


I 






Uii______---^; 

^'//^ 






167rqa>l 


//VI 


/ VII 




167rqa=l 


c=o / 

(V) / 


C<0 




r 

167rqa<l 





FIG. 4: Qualitative behaviour of A(r) for the IRD-A2 models 
when a > 1 . This behaviour is similar in the three cases b), 
c), a), {16-nqa = 1, respectively), excepting by the sign of the 
metric function of the critical configurations at the center, 
which is \ C rit(r = 0) = 0, respectively. Then, by simplicity, 
we have plotted a unique set of curves for the different con- 
figurations and modified the relative position of the A = 
axis for each case (continuous horizontal lines). The metric 
functions diverge at the center to +oo when C > and to 
— oo for C < 0, whereas for the critical configuration (C = 0), 
they attain a finite value, which is positive, null or negative, 
depending on the case. The slope of the critical configuration 
at the center diverges for this value of a. Curve I corresponds 
to naked singularities in all cases. Curve II corresponds to 
extreme black hole configurations in the 16nqa > 1 case and 
to naked singularities otherwise. Curves III and IV corre- 
spond to two-horizons black holes in the 16ivqa > 1 case and 
to naked singularities otherwise. The dashed curve V is the 
critical configuration and corresponds to single-horizon black 
holes in the 16nqa > 1 case, to extreme black points when 
16nqa = 1 and to naked singularities in the 16nqa < 1 case. 
Curves VI and VII correspond to single-horizon black holes in 
all cases. Curves IV and VI show how the metric approaches 
the behaviour of the critical case as C — > from above and 
below, respectively. 



<p(X,Y = 0,n) 



In 



1 — [i\/X 



- 2VX 



for the restriction for purely electric fields of the family 
of Lagrangian densities supporting ESS solutions of the 
form (141 [) . This expression is valid in the range < 
X < l//i 2 , where it behaves as expected for these 1RD- 
A2 cases. The energy density for these solutions reads 
now as 



p{r, Q,n) = 



2a/I + R 2 , / Vl + R 2 + 1 
— In 



R 2 



VI + R 2 - 1 



(44) 
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A(r) 




general statements a) to c) for the different subcases of 
this example are straightforward. 



C. IRD-UVD case 

In this case the internal and external integrals of energy 
do not exist, but the primitive of r 2 T$ (r, Q) in equation 
(1261) is well defined up to an arbitrary constant 

s(r, Q, r) = 4tt J r 2 T °(r, Q) + Y. (46) 
Thus the metric is integrated as 



FIG. 5: Same comments as in Fig. 4, but now for a < 1 . The 
slope at the center of the metric function for the critical con- 
figuration (C = 0) vanishes in this case (see also the discussion 
in the text). 



A(r) 




FIG. 6: Same comments as in Fig. 4, but now for a = 1 . The 
slope at the center of the metric function for the critical con- 
figuration (C = 0) is finite in this case (see also the discussion 
in the text). 



\{r,Q,D) = 1 + - 
r 



D 2e(r,Q,0) 



(47) 



where we have introduced the arbitrary constant D = 
C — 2T. Owing to the admissibility conditions and the 
previous analysis, the function r 2 TQ(r, Q) is a monotonic 
function, decreasing from +oo at r = 0, where it behaves 
as r p (p < —I), and vanishing asymptotically as r q (— I < 
q < 0) (see Eqs. flTTj) ) and (|T2"j) ). Consequently, the curves 
e(r, Q, r) are monotonic and convex, increasing from — oo 
at r = and exhibiting a horizontal parabolic branch as 
r — > oo (for any V and, in particular, for F = 0). The 
condition for the vanishing of the metric function is now 



e(r,Q, r = 0) 



D 



(48) 



whose solutions define the radii of the horizons of the 
different black hole configurations (see Fig. 2). There is 
always a value of the constant D in the r.h.s. of Eq. (|4"8")l 
determining the straight line of the beam ^^p- which is 
tangent to the energy function and defining an extreme 
black hole, whose associated radius ri iextr {Q) is the solu- 
tion of 



where R = — -, — 
function becomes 



for any sign of Q. The internal-energy 



e in {r,Q,n) 



8^Q 3 / 2 



Ry/l + R 2 



-R A In 



Vi + R 2 + i 



VT+W-i 



+2\n(R+ ^Jl + R 2 ) , 



(45) 



which exhibits the qualitative behaviour of the A2-curves 
of Fig. 2. From these formulae the determination of the 
structure of the metrics and the verification of all the 



47rr 2 T °(r,Q) 



(49) 



(note that, as in the previous cases, this radius does not 
depend on the value of the constant D = C — 2T). Equa- 
tion (|48| allows for the determination of a value of the 
constant D extr (Q) corresponding to this tangent line. We 
see that, aside from extreme black holes, there are 
now configurations with a naked singularity at the cen- 
ter (when D > D extr (Q)) and also two-horizons black 
holes (when D < D extr (Q))- Excepting by the anoma- 
lous asymptotic behaviour, these structures are qualita- 
tively similar to those associated to the RN solutions. 

The dominating term in the expression of the metric 
function near the center and as r — ¥ oo is the last one 
of Eq.(|4"Tf and, consequently, this function diverges to 
+oo at r = and behaves asymptotically as I — A(r —¥ 



10 



oo, Q, D) ~ r q —> + (no matter the sign of the constant 
D), exhibiting a similar form as that of the curves I to 
IV in Fig. 3, with a single minimum \{r min ,Q,D) < 1. 
As can be easily checked by deriving Eq. (jT71) and taking 
into account the admissibility conditions, the values of 
these minima decrease without limit and their locations 
approach the origin as D — > — oo. 

As an illustrative example of this class of models let us 
consider the field 



r 



(50) 



which is the ESS solution of a very simple IRD-UVD 
model, but exhibiting all the features of this case. The 
restriction to the half axis X > 0, Y = of the corre- 
sponding family of Lagrangian densities supporting this 
solution is obtained from Eq.© and takes the form 



<p(X,Y = 0) = -A 3 / 2 . 



The energy density reads 



4Q 



3/2 



3r 3 ' 



and the primitive of r 2 Tg (r, Q) is 

, 16ttQ 3 / 2 
£[r,Q) = g ln(r) 

The metric function takes the form 



(51) 



(52) 



(53) 



(54) 



w D 32ttQ 3 / 2 ln(r) 

X(r,D,Q =1 + J 

r 3 r 

which diverges to +00 at r = 0, reaches a minimum with 
A < 1 and approaches asymptotically to A — > l - as r — >• 
00, slower than the Schwarzschild solution. The extreme 
black hole horizon radius is 



rhextr(Q) = 



32irQ 3 / 2 



(55) 



and the corresponding value of the constant D extr {Q) is 



Dextr{Q) — 



32ttQ 3 / 2 



.»(^i-, 



(56) 



As easily verified, the solutions (|5"1|) with D > D extr (Q) 
have minima X m in > and are naked singularities. Those 
with D < D extr {Q) have X m in < 0, cutting twofold the 
A = axis, and correspond to two-horizons black hole 
configurations (see Fig. 3, curves I to IV). 

The lagrangian ([ST]) and the associated solutions ([50j) 
correspond, respectively, to the limit case of the IRD-Al 
family (|43|) and their associated solutions (|4Tj) as the pa- 
rameter fi —} 0. In this limit the maximum field strength 
(a = /i _1 ) of the ESS solutions of this family diverges 
and the limit model becomes IRD-UVD. 



IV. SUMMARY AND PERSPECTIVES 



Let us summarize the main results obtained until now. 

The contents of this paper, together with those of 
Ref.[l2j], complete the analysis of the geometrical struc- 
ture of the G-ESS solutions of admissible NEDs min- 
imally coupled to gravity. The families of Lagrangian 
theories considered exhaust the class of the admissible 
ones. 

These admissible G-NEDs lead always to asymptoti- 
cally flat G-ESS solutions. At large distances, the met- 
rics approach flatness as the Schwarzschild field in cases 
Bl, B2 and B3, with a well defined ADM mass. On 
the other hand, they approach flatness slower than the 
Schwarzschild field for models with IRD asymptotic be- 
haviour and the ADM mass is not defined in these cases. 
The central-field behaviour plays no role in determining 
the asymptotic gravitational structure of the solutions. 

The geometrical structures of the G-ESS solutions at 
finite r are dependent on the central field behaviour of 
the ESS solutions in flat space or, equivalently, on the 
behaviour of the Lagrangian densities on the boundary 
of their domain of definition at Y = 0. The present 
analysis, together with the one of Ref . [I2I] . show that 
the qualitative nature of these structures is indepen- 
dent on the asymptotic behaviours of the G-ESS solu- 
tions. For the UVD central field cases, combined with 
Bl, B2 or B3 behaviours as r —> 00, these structures 
are similar as those of the Reissner-Nordstrom solutions 
of the Einstein-Maxwell field equations (naked singular- 
ities, extreme black holes and two- horizons black holes). 
The UVD-IRD cases exhibit similar structures, but with 
an anomalous (non Schwarzschild-like) asymptotic be- 
haviour. For Al and A2 central field cases there are in ad- 
dition single-horizon black holes and (in case A2) extreme 
and non-extreme black points, no matter the asymptotic 
properties (and, hence, the finite or divergent- energy 
character) of the ESS solutions in flat space. But when 
combined with IRD behaviour at large r, they lead also 
to anomalous asymptotic behaviours for the gravitational 
fields. 

As already mentioned, all the G-ESS solutions con- 
sidered here have a curvature singularity at the center, 
whose causal structure is dependent on the sign of the 
metric function there. For NEDs with Maxwellian weak- 
field limit (B2 cases) this is a consequence of the Bron- 
nikov theorem [l3| . It is easy to extend the validity of this 
theorem for the admissible models belonging to the other 
asymptotic cases (Bl, B3 and IRD) considered here. Let 
us summarize the argument of Bronnikov, as presented 
in Ref. [13: 

Owing to the diagonal character of the components 
of the Ricci tensor in the coordinate system (l24l) . the 
quadratic curvature invariant RapR 01 ® can be written as 
the sum of squares of the diagonal components of the 
mixed Ricci tensor 
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(57) 



As established in reference [TT| the stability of the ESS 
solutions of these models in flat space is endorsed by the 
necessary and sufficient condition 



Then, the finiteness of this invariant at a given point 
implies the finiteness of the diagonal components of the 
Ricci tensor there and, consequently, of those of the 
energy-momentum tensor (|23l) . in such a way that the 
quantity 



Y=0 



(58) 



must be finite. On the other hand, the first-integral 
gives 



E 



dtp 
dX 



Q 



Y=0 



(59) 



At the center, the r.h.s. of (|59| diverges, and the compat- 
ibility of this expression with the finiteness of ([55)) leads 
to the conditions 



dip 
dX 



X^0,Y=0 



E(r -> 0) -> 0. 



(60) 



This behaviour is incompatible with models belonging to 
the Maxwellian cases B2, which proves the theorem in its 
initial version. 

The behaviour (f6"0"|) is also incompatible with the cases 
IRD and Bl, for which the theorem is automatically ex- 
tended. The first one of conditions (|60| is compatible 
with the cases B3 and, at first sight, one could expect the 
existence of regular G-ESS solutions for some G-NEDs of 
this class. However, the positivity of the energy requires 
for the ESS solutions that E(r — >• oo) —> 0. This con- 
dition, aside from the second one of Eq. (|6"0")) (and the 
positivity of ^§ |y_ ) breaks the monotonic character of 
E(r). As a consequence, the first-integral (1591 implies 
the multi-branched character of tp(X, Y = 0) and the 
non-admissibility. This completes the extension of the 
Bronnikov theorem: The G-ESS solutions of any admis- 
sible NED minimally coupled to gravity exhibit a curva- 
ture singularity at the center. 



dx - dv^ 



(61) 



to be satisfied by the Lagrangian densities in the domain 
of definition of the solution. The question of whether 
ESS solutions of admissible models satisfying this sta- 
bility flat-space criterion are also stable when minimally 
coupled to gravity or, on the contrary, the gravitational 
coupling has the effect of destabilizing stable configura- 
tions in flat space (requiring new conditions to be im- 
posed for stability in the gravitational case) is currently 
under analysis. 

In this work we have not addressed the problem con- 
cerning the thermodynamic analogies for these models. 
As it has been established in Ref. [19J the zeroth and first 
laws of black hole thermodynamics can be extended to 
any G-NED with Maxwellian weak-field limit (B2 cases) 
(although the Smarr formula cannot, in general). It is 
easy to extend this result to cases Bl and B3 and it is 
possible to obtain generalized Smarr formulae for admis- 
sible G-NEDs. These and other related issues will be 
analyzed elsewhere. 

Finally, although it lies beyond the scope of this pa- 
per, let us mention, as another important issue in this 
context, the study of the dynamics of solutions and wave 
propagation in G-NEDs. Among the first analysis involv- 
ing this kind of theories, the one of Ref.[2(| generalizes 
to the gravitating case, for Bl electrodynamics, the Boil- 
lat results on wave propagation of general NEDs in flat 
space 21] . More recent advances in this line of work 
with applications to modern problems can be found, for 
example, in Ref. (22| . 
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